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FINAL  REPORT:  Grant  #  W911NF-13-1-0474 


Numerical  simulation  of  atmospheric  boundary  layer  flow  over 
battlefield-scale  complex  terrain:  surface  fluxes  from  resolved  and  subgrid 

scales 

A  final  report  submitted  to  the  Army  Research  Office,  Environmental  Sciences 
Directorate,  Atmospheric  Sciences  Program  (Grant  #  W911NF-13-l-0f7f;  PM:  Dr.  S. 

Collier,  Dr.  G.  Videen) 

PI:  William  Anderson,  Mechanical  Engineering  Department,  University  of  Texas  at  Dallas 

1  Background 

This  document  presents  a  summary  of  activities  during  Year  1  of  Grant  ff  W911NF- 
13-1-0474  (FY2014).  The  funding  period  began  in  September,  2013,  when  the  PI  was  at 
Baylor  University.  As  discussed  more  extensively  in  Section  6,  the  PI  moved  from  Baylor 
to  the  University  of  Texas  at  Dallas  (UTD)  in  July,  2014.  This  professional  move  has 
affected  the  project  by  forcing  the  “Basic  Award”  to  in  fact  conclude  at  Baylor,  while  a 
new  start  basic  award  will  begin  at  UTD  in  October,  2014,  for  a  duration  of  two  years. 
This  notwithstanding,  the  Year  fi  1  research  effort  has  been  remarkably  successful.  When 
the  original  proposal  was  submitted  to  then-PM  Dr.  Gorden  Videen  (Spring,  2012),  PI 
Anderson  had  proposed  to  develop  numerical  treatments  for  fluxes  of  momentum  and 
scalars  in  turbulent  atmospheric  surface  layer  flows  over  complex  topographies  at  “bat¬ 
tlefield  scales”  ~  (9(10  km).  This  work  was  effectively  following  on  from  Anderson’s  PhD 
studies  at  Johns  Hopkins  University  [1,  2],  The  work  was  motivated  by  limitations  of 
present  computing  resources  for  modeling  microscale  atmospheric  flows  (with  large-eddy 
simulation,  LES)  and  accounting  for  the  full  spectrum  of  topographic  modes  typical  of 
a  natural  terrain.  For  example,  Figure  1(a)  shows  a  digital  elevation  map  from  Mason 
County,  Texas,  retrieved  from  the  United  States  Geological  Survey  seamless  server.  The 
horizontal  domains  are  ( Lx,Ly )  «  (10  km,  10  km).  The  fine-grained  height  field  (Figure 
la)  is  resolved  at  Af  =  3  m.  LES  with  resolution  sufficiently  fine  to  represent  these 
small  topographic  scales  while  still  capturing  the  largest  ABL  turbulent  structures  would 
require  significant  (perhaps  prohibitively  large)  scientific  computing  resources.  Consider 
then  the  effect  of  filtering  the  height  field  at  a  scale  more  appropriate  in  LES  with  this 
domain  size,  Ales  —  96  m,  Figure  lb.  These  images  demonstrate  that:  ( i )  filtering 
removes  important  information  from  the  height  field;  (ii)  the  landscape  roughness  is 
heterogeneous;  and  (Hi)  fluvial  landscapes  are  characterized  by  branching,  fractal-like 
channel  networks  [3].  We  note  that  the  choice  of  location  for  Figure  1  is  arbitrary.  While 
Anderson’s  preliminary  efforts  [1,  2]  have  focused  on  complex  topographies  most  like  that 
shown  in  Figure  1,  one  could  also  imagine  addressing  for  example  flows  over  partially  re¬ 
solved  urban  environments [4,  5,  6],  or  vegetative  canopies[7,  8,  9,  10]. 

As  a  PhD  student  at  Johns  Hopkins,  Anderson  worked  with  PhD  adviser  (C.  Meneveau) 
to  develop  models  for  momentum  fluxes  (that  is,  aerodynamic  drag)  associated  with 
unresolved  topographic  modes  in  the  Figure  1(b)  topography.  This  was  accomplished  by 
first  supposing  that  the  total  aerodynamic  drag  (i.e.  due  to  h(x,y ))  could  be  decomposed 
into  that  imposed  by  the  filtered  ( h(x,y ))  and  unresolved  (h'(x,y))  components: 

FA  =  ~  [  Pw  hi  dS  +  p  f  t%a  hj  dS ,  (1) 

Js  Js 
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where  ht  is  the  unit  vector  normal  to  the  hltered  landscape,  h  (x,  y ),  and  pw  is  the  resolved 
wall  pressure  acting  on  h(x,y),  while  drag  forces  due  to  h'  ( x,y )  are  represented  by  the 
equilibrium  logarithmic  law: 

w.  A  ~  r  tt 

R3  2^  =  _  KlJ 
P  UtU  [\og(z/z0) 

where  Zo  is  a  momentum  roughness  length,  n  is  the  von  Karman  constant,  uT  is  the  shear 
velocity,  and  U  is  magnitude  of  the  horizontal  components  of  the  velocity  vector*.  Next, 
we  invoke  a  self-consistency  argument  (the  Gernrano  identity [11])  on  the  plane-averaged 
aerodynamic  drag  at  a  grid-  and  test-filter  scale;  this  is,  we  enforce  that  the  total  drag 
embodied  in  Equation  1  should  be  invariant  to  the  filter  resolution  (assuming  that  the 
filter  range  is  in  the  “inertial”  range  of  the  topographyfl,  2]).  Thus,  one  could  imagine  a 
test-filtered  version  of  Equation  1: 

if A  =  -  [  v  hi  dS  +  p  f  t%’2A  hj  dS ,  (3) 

Js  Js 

where  bG  denotes  a  test-filtered  (more  smooth)  quantity.  The  self-consistency  argument: 

(Fa)  =  (if  A>,  (4) 

states  that  the  averaged  aerodynamic  drag  is  equivalent,  regardless  of  the  choice  of  filter 
width,  A.  Thus,  as  A  is  varied,  the  drag  imposed  by  the  subgrid  and  resolved  scales  should 
change  systematically  such  that  the  net  quantity  is  preserved.  Equation  4  is  insightful, 
since  it  also  allows  one  to  inform  parameters  needed  to  parameterize  the  subgrid-scale 
roughness  length,  z0.  For  brevity,  additional  details  regarding  the  Equations  1  to  4  de¬ 
velopment  are  excluded  here  (Anderson  has  presented  this  work  many  times,  including 
during  an  invited  talk  at  the  Army  Research  Laboratory  Headquarters  in  Adelphi,  MD 
in  July,  2013). 


u 

U' 
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It  was  originally  anticipated  that  the  1  to  4  development  would  be  applied  analogously 
to  model  surface  fluxes  of  passive  scalars  for  Figure  1-like  complex  topographies.  As 
a  first  stage  effort,  we  sought  to  study  the  passive  scalar  surface  flux  characteristics 
of  flow  over  a  fully-resolved  fractal-like  topography.  The  equilibrium  logarithmic  law 
expression  for  passive  scalar  fluxes,  q"  (neutral  stratification  -  stability  correction  terms 
not  needed[12,  13])  is: 


=  9tut 

pCp 


K  (9s  -  0) 
\og(z/z0S)UT' 


(5) 


where  9T  is  shear  temperature,  20s  is  scalar  roughness  length,  cv  is  specific  heat,  9S  is  the 
surface  temperature,  and  9  is  passive  quantity  (for  example,  temperature  in  the  absence 
of  large  thermal  gradients).  A  large  body  of  literature  exists  regarding  values  of  z0s 
[14,  15,  16,  17,  18,  19,  20,  21]  in  terms  of  the  momentum  roughness  length, 


z0s  =  zo  exp 


— k  (stp 1  -  1/2)  , 


(6) 


*In  this  report,  we  adopt  the  following  nomenclature:  x  =  {x,y,z}  correspond  with  the  streamwise, 
spanwise,  and  vertical  positions,  respectively;  u  =  {u,v,  w}  correspond  with  the  streamwise,  spanwise, 
and  vertical  velocity  components,  respectively. 
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Figure  1:  Contour  view  of  height  distribution  at  30.57°  :  30.77°  N  and  —99.22°  : 
—99.47°  E  (approximately  Mason  County,  Texas).  Part  (a)  illustrates  unfiltered  height 
held,  h  (x,y),  resolved  at  l/9th  Arc-Seconds.  Part  (b)  illustrates  filtered  height  held  after 
filtering  at  LES  computational  mesh  resolution,  h(x,y). 


where  Sto  and  Cd$  are  the  interfacial  Stanton  number  and  interfacial  drag  coefficient. 
The  term  within  the  exponential  in  Equation  6  is  often  lumped  into  the  parameter, 
B -1  =  St^1  —  CcIq1^  [22,  23],  and  named  the  interfacial  transfer  coefficient.  Numerous 
models  for  B have  been  developed,  all  of  which  have  a  dependence  on  the  roughness 
Reynolds  number,  Re0  =  uTz0/u,  as:  St^1  =  (ftiRep  Pr^3 .  Various  groups  have  offered 
values  for  the  set  of  (j)  parameters,  while  Pr  is  the  Prandtl  number.  In  a  recent  study, 
we[24]  used  an  immersed  boundary  method[25]  to  model  how  over  fully-resolved,  fractal 
topographies.  Since  the  topographies  were  fully-resolved,  we  could  dynamically  evaluate 
the  equivalent  zo  during  simulation,  thereby  allowing  a  priori  computation  of  Re o,  Sto, 
z0s  (Equation  6),  and  then  q"  (Equation  5).  The  principle  finding  from  this  study [24]  was 
that  the  ratio,  varied  over  several  orders  of  magnitude  as  Re0  changed;  moreover, 

z0 s/ Zo  «  1CT1  for  most  of  the  fractal  topographies  considered.  This  conflicts  somewhat 
with  the  commonly-used  assumption  that  ^os/V)  =  10_1,  and  is  consequential  for  numer¬ 
ical  weather  prediction  of  micrometeorological  conditions  over  battlefield  environments. 
This  work  was  largely  completed  during  the  interval  between  original  proposal  submission 
and  eventual  award  (Spring,  2012  to  Summer,  2013).  By  the  time  the  grant  was  awarded, 
Anderson  had  developed  a  strong  interest  in  studying  the  morphological  characteristics  of 
turbulent  atmospheric  flows  over  urban  environments.  At  the  PM’s  approval,  Anderson 
effectively  allocated  Year  ^  1  of  the  funding  effort  to  this  problem,  outlined  below. 

2  Findings 

Turbulent  momentum  transport  in  atmospheric  boundary  layer  flows  is  of  critical  im¬ 
portance  to  the  performance  of  modern  wind  farms [26],  aerodynamics  of  vegetative 
canopies[27,  9]  and  urban  environments [28,  29,  5,  4],  and  geomorphological  processes 
associated  with  evolution  of  aeolian  desert  landscapes [30].  These  examples  are  distinctly 
different  to  smooth  wall  turbulence,  owing  to  the  presence  of  a  distribution  of  obstacles 
of  height,  h,  that  protrude  into  the  inertial  layer  of  the  flow.  These  obstacles  absorb 
momentum  through  pressure  drag,  induce  flow  separation  and  serve  to  produce  obstacle- 
scale  coherent  motions  that  occupy  the  region  between  the  wall  and  two  to  four  times  the 
obstacle  height  -  the  roughness  sublayer[31,  6].  In  the  roughness  sublayer,  the  vertical 
aerodynamic  drag  distribution  due  to  the  presence  of  obstacles  results  in  an  inflected 
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mean  streamwise  velocity  profile  at  the  approximate  average  obstacle  height  [8,  32],  As  a 
result,  the  mean  flow  gradient  exhibits  its  maxima  at  the  inflection [30]  (not  the  wall),  and 
the  turbulence  kinematics  are  fundamentally  different.  For  flow  over  vegetative  canopies, 
Raupach  et  al.  [33]  illustrated  that  flow  in  the  roughness  sublayer  resembled  a  turbulent 
mixing  layer  with  positively  and  negatively  skewed  streamwise  and  vertical  velocity  fluc¬ 
tuations,  respectively,  and  Reynolds  stresses  composed  predominately  of  “sweeps” .  With 
this,  the  turbulence  morphology  is  characterized  by  Kelvin-Hclmholtz  spanwise  “rollers” 
which  originate  at  the  velocity  profile  inflection  and  undergo  a  downstream  metamorpho¬ 
sis  leading  to  hairpin  packets  [27].  Recently,  Ghisalberti[34]  demonstrated  the  existence 
of  an  universality  in  roughness  sublayer  statistics  for  flows  over  diverse  canopies,  and  in¬ 
troduced  the  term  “ obstructed  shear  flow ”  to  categorize  such  flows.  Within  the  sublayer, 
the  turbulence  geometric  macroscale  is  effectively  set  by  the  vorticity  thickness[35,  33], 
8W  =  Uh/(dU/dz)\h,  and  streamwise  spacing  of  vortex  cores,  Ax  =  2nL1Xl  where  L is 
the  integral  length  scale.  The  mixing  efficiency  associated  with  these  Kelvin-Hclmholtz 
roughness  sublayer  motions  (sustained  by  mean  flow  gradient  with  maxima  at  the  canopy 
height)  exceeds  the  analogous  value  for  a  logarithmic  profile  by  approximately  50  %. 
Thus,  turbulence  in  the  roughness  sublayer  is  characterized  by  vigorous  mixing  and  com¬ 
plex  structural  attributes. 

Above  the  roughness  sublayer  -  in  the  inertial  layer  -  the  mean  flow  profile  exhibits  log¬ 
arithmic  scaling  and  the  turbulence  structure  resembles  smooth  wall  flow  [8,  32,  5].  That 
is,  the  domain  is  occupied  by  persistent  streamwise-elongated  coherent  parcels  of  rela¬ 
tively  low  and  high  momentum  regions.  For  smooth  walls,  such  low  momentum  regions 
meander[36],  and  are  encapsulated  by  hairpin  packets[37,  38,  39,  40,  41,  42,  43,  44]  at 
the  interface  between  zones  of  quasi-uniform  momentum[45].  For  rough  walls,  the  pres¬ 
ence  of  persistent  structures  is  also  reported [5,  42],  however  experiments  have  indicated 
that  h- scale  coherent  motions  associated  with  the  roughness  sublayer  seemingly  attenu¬ 
ate  the  lengths  of  logarithmic-layer  coherent  motions[42].  Recently,  Coceal  et  al.[ 5]  used 
direct  numerical  simulation  (DNS)  to  study  flow  over  a  staggered  array  of  uniform  height 
cubes  [28]  with  characteristic  scale,  h/H  —  where  H  is  the  boundary  layer  depth.  They 
illustrated  the  existence  of  hairpin  packets  (and  “cane”  structures:  inclined  coherent  par¬ 
cel  with  only  one  leg  of  the  hairpin[41]  around  the  low  momentum  region).  Moreover, 
they  used  approximated  conditional  averages [46,  37]  to  illustrate  the  significant  reduc¬ 
tions  in  streamwise  coherence  of  their  rough  wall  turbulence  relative  to  a  smooth  wall. 
These  efforts  to  characterize  the  structural  attributes  of  smooth  and  rough  wall  turbulence 
have  typically  focused  on  spatial  characteristics.  Further,  most  previous  studies  address 
either  instantaneous,  ensemble-averaged,  or  time-averaged  statistically  stationary  turbu¬ 
lence  statistics.  Here,  we  sought  to  reconcile  spatial  and  temporal  turbulence  statistics 
during  flow  over  urban-like  topographies,  and  topographies  with  a  uniform  aerodynamic 
roughness  length.  This  was  accomplished  simply  by  nominating  locations  throughout  the 
computational  domain  at  which  flow  statistics  are  recorded  across  the  depth  of  the  bound¬ 
ary  layer  in  time  (thus,  we  use  LES  to  generate  datasets  that  would  more  typically  be 
associated  with  held  campaign  data  from  a  tower  equipped  with  sonic  anemometers)  [47] . 

We  used  large-eddy  simulation  (LES)  to  model  how  over  topographies  composed  of  a 
staggered  distribution  of  wall-mounted  cubes  (one  of  the  cases  considered  is  “C20S”  from 
the  wind  tunnel  study  by  Cheng  and  Castro [28],  and  considered  in  the  more  recent  DNS 
study  by  Coceal  et  al.  [5] ) .  Figure  2(a)  shows  a  perspective  view  of  the  array  of  cubes 


4 


U0 


Figure  2:  Illustration  of  [28]  block  case  (“C20S”)  considered  as  lower  topography  in 
present  study:  (a)  perspective  image  showing  freestream  flow  direction,  Uq\  and  (b)  plan- 
form  image  with  indication  of  positions  at  which  time-height  velocity  data  are  recorded. 
Owing  to  the  topography  attributes,  the  four  points  effectively  capture  transient  dynam¬ 
ics  at  all  points  in  the  domain.  The  above  topography  is  labeled  here  Case  SC2. 
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Figure  3:  (Color)  Illustration  of  computational  domains  considered  in  the  present  study: 
(a)  uniform  height,  staggered  array  of  blocks[28]  with  h  varying,  s/h  =  1;  and  (b)  ho¬ 
mogeneous  roughness,  zq.  For  Panel  (a),  see  also  Figure  2.  For  all  cases,  H  =  1000 
m.  Illustrative  locations  of  positions  at  which  time-series  of  streamwise  and  vertical  ve¬ 
locity  are  recorded  are  shown  for  discussion  (vertical  red  lines).  In  both,  an  illustrative 
time-averaged  streamwise  velocity  profile  is  shown,  where  Uq  is  the  “outer”  freestream 
velocity. 


with  edge  length,  h.  Figure  2(b)  shows  the  arrangement  in  planform,  with  indication  of 
the  streamwise  spacing,  s/h,  between  rows  of  cubes  (here,  we  consider  only  cases  with 
s/h  =  1).  Figure  2(b)  also  includes  indication  of  Points  X\,  x2,  x3,  and  x 4;  during  LES, 
time-series  of  streamwise  and  vertical  velocity  have  been  recorded  across  the  depth  of  the 
boundary  layer  at  these  points.  This  is  better  illustrated  in  Figure  3(a),  which  shows  the 
cubes  and  streamwise  spacing  (for  the  cube  cases,  the  y  coordinates  of  Positions  X\  to  X4 
are  equal  and  set  to  intersect  the  center  of  cubes;  the  x  coordinates  are  varied  such  that: 
X\  is  precisely  at  the  center  of  a  cube,  x2  is  centered  between  the  upwind  and  downwind 
row  of  cubes,  x3  is  precisely  at  the  center  of  adjacent  cubes,  and  x4  is  centered  between 
the  upwind  and  downwind  row  of  cubes).  For  generality,  we  also  model  flow  over  homo¬ 
geneous  roughness  lengths,  shown  in  Figure  3(b).  This  work  is  novel,  since  it  provides 
datasets  that  might  more  readily  be  retrieved  from  a  micrometeorological  tower  equipped 
with  sonic  anemometers  (while  retaining  the  spatial  flexibility  offered  from  LES).  From 
adopting  such  a  viewpoint,  we  can  gain  deeper  insight  about  the  role  of  coherent  struc¬ 
tures  in  the  observed  time-series  phenomena. 

Figure  4  shows  time- height  contour  maps  of  fluctuating  velocity  component  (defined  for 
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this  purpose  as  u'(x,t)  =  u(x,t)  —  ( u(x,t))t ,  where  (. .  ,)t  denotes  averaging  over  time, 
t).  For  brevity,  we  show  only  u'  and  w'  at  the  limiting  positions  (i.e.  X\  and  £4),  although 
intermediate  values  exhibit  similar  patterns.  We  also  show  the  same  contours  for  flow 
over  a  homogeneous  roughness  distribution  (Figures  4e  and  4f).  It  is  clear  from  inspection 
that  the  u'  contours  exhibit  a  significant  advective  lag  for  all  cases;  that  is,  the  passage  of 
parcels  of  relatively  high-  or  low-momentum  in  the  aloft  boundary  layer  precedes  “simi¬ 
lar”  activity  in  the  roughness  sublayer.  This  is  further  evidenced  by  consideration  of  the 
solid  contours,  which  illustrate  Reynolds  stresses  owing  to  “sweep”  and  “ejection”  events 
(determined  via  quadrant  analysis [48]).  One  can  further  appreciate  that  -  in  general  - 
elevated  Reynolds  stresses  owing  to  sweeps  (vigorous  downward  excursions  of  high  mo¬ 
mentum  fluid,  into  the  sublayer)  are  also  preceded  by  aloft  passage  of  high  momentum 
regions.  This  is  clear  evidence  that  the  characteristics  of  outer  layer  coherent  parcels  of 
fluid  are  “imprinted”  on  roughness  sublayer  processes;  moreover,  qualitative  observations 
from  Figures  4(a),  (c),  and  (e)  suggest  an  advective  lag  between  passage  of  such  outer 
layer  structures  and  excitation  of  sublayer  dynamics. 

We  quantify  the  u'  advective  lag  qualitatively  observed  in  Figure  4.  This  is  accomplished 
in  two  steps.  Firstly,  a  reference  height,  ZRef.,  is  nominated  at  which  a  reference  dataset 
can  be  collected.  We  selected  znet\  to  be  slightly  above  the  “top”  of  the  roughness. 
We  adopted  this  approach  since  it  facilitates  comparison  between  these  fundamentally 
different  topographies.  From  here,  we  compute  correlations  maps  as  the  convolutions: 

/OO 

u  (xi,z,t)u(xi,zRe{.,t  +  r)dt ,  (7) 

-OO 

where  Xi  represents  a  spatial  position  (Figure  2b).  From  this,  we  compute  advective  lag 
as: 

Tmax.(z)  =  argmax(('U/(a3;,  Z,  t)  *  u'(Xi,  ZRef.,  t))(t))i  (8) 

t. 

where  we  discard  rmax.  values  corresponding  with  7 (z,rmax.)  <  X >  where  y  —  0-3  is  a 
predefined  threshold.  We  experimented  with  a  range  of  y  values,  hireling  generally  that 
increasing  the  threshold  only  serves  to  remove  spurious  values  for  z  »  ZRef.  while  the 
underlying  rmax.(^)  trends  were  robust  and  indifferent  to  y.  Figure  5  shows  shear  nor¬ 
malized  advective  lag,  Tmax.(z)ur.ff’_1  for  Cases  SCI,  SC2  and  U1  (Cases  SCI  and  SC2 
correspond  with  h/H  =  1/8  and  =  1/4,  respectively).  For  Cases  SCI  and  SC2,  it  is  clear 
that  Positions  X\  to  £4  exhibit  effectively  the  same  rma.T  ii-r II  1  profiles.  Above  ZRef. ,  the 
profiles  are  roughly  linear,  and  we  remind  the  reader  that  2Ref.  is  the  approximate  obsta¬ 
cle  elevation  (or  center  of  the  mean  streamwise  velocity  profile  inflection).  Further,  the 
advective  lag  is  always  negative  above  ZRef.,  which  is  precisely  consistent  with  qualitative 
observations  in  Figure  4.  For  Position  aq,  there  of  course  are  no  Tmax.rtTi7_1  datapoints 
below  ,2 Ref.  due  to  the  solid  cube.  A  few  datapoints  are  available  for  Position  37 ,  owing  to 
its  position  in  the  lee  of  the  cubes.  However,  at  Positions  £3  and  £4,  there  are  rmax  ii^II  1 
values  which  illustrate  that  activity  at  2Ref,  precedes  activity  below  h  and  is  positively 
correlated  with  momentum  transport  processes  in  the  canopy,  not  the  roughness  sub¬ 
layer.  We  emphasize  also  that  the  TmaxuTH~l  values  in  Figure  5  are  only  associated 
with  7  exceeding  threshold,  y  =  0.3.  Thus,  the  approximate  linearity,  rmax.  ~  —  z  for 
z  >  2Ref.,  is  a  product  of  actual  processes  within  the  roughness  sublayer  and  inertial  layer. 

The  presence  of  meandering,  coherent  parcels  of  relatively  low  and  high  momentum  in 
turbulent  wall-bounded  flows  over  smooth[37,  38,  39,  40,  41,  43,  44]  and  rough[5,  42]  walls 
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is  well  known.  The  low  momentum  regions  (LMR)  are  encapsulated  by  hairpin  packets 
at  the  interface  between  regions  of  differing  momentum [45,  41].  For  the  case  of  cube 
roughness  such  as  the  cases  considered  here,  Figure  6  is  a  sketch  of  the  aforementioned 
structural  attributes.  A  streamwise-vertical  transect  through  a  LMR  would  reveal  a  typi¬ 
cal  inclination  of  7  ^  15°  [37,  42],  Quantitative  visualization  of  the  Figure  6  dynamics  can 
also  be  attained  with  vortex  visualization  techniques  (excluded  here,  for  brevity).  Figure 
5  shows  vertical  profiles  of  advective  lag,  Tmax  urH~l,  based  on  a  reference  elevation  just 
above  the  canopy  (see  Figure  5).  The  linearity  exhibited  by  rmaY  u^H1  points  to  an 
underlying  physical  process  in  the  roughness  sublayer  and  inertial  layer.  We  attribute 
Tmax.  ~  —  z  scaling  to  the  passage  of  regions  of  alternating  low  and  high  momentum  in 
the  roughness  sublayer  and  inertial  layer. 

Following  the  passage  of  a  LMR,  the 
first  position  in  the  domain  to  ex¬ 
perience  relatively  higher  momentum 
would  be  in  the  inertial  layer,  above 
the  LMR.  As  the  LMR  advects  down¬ 
stream,  relatively  higher  momentum 
would  be  recorded  at  progressively 
lower  elevations.  I11  keeping  with 
the  viewpoint  adopted  thus  far,  time- 
series  measurement  of  flow  statistics 
throughout  the  domain  at  a  position 
upwind  of  the  LMR,  27,  would  experi¬ 
ence  low  momentum  higher  in  the  do¬ 
main  first.  The  LMR  would  be  ob¬ 
served  at  progressively  lower  z  with 
increasing  time.  The  above  pattern 
would  be  true  for  regions  of  relatively 
high  momentum  too.  Since  represen¬ 
tative  information  is  known  about  the 
macroscale  attributes  of  these  coher¬ 
ent  motions,  we  have  developed  a  sim¬ 
ple,  semi-empirical  model  to  predict 
the  advective  lag  between  passage  of  such  motions  in  the  inertial  layer  and  evidence 
of  their  “imprint”  on  roughness  sublayer  and  canopy  dynamics. 

If  a  LMR  (i.e.  u'  <  0,  above)  has  depth,  5',  its  length  can  be  evaluated  based  on 
assumption  of  a  hairpin  train  inclination  angle: 

Ls  «  5'/  tan  (7) .  (9) 

Furthermore,  if  we  assume  that  a  representative  advective  velocity  for  low  and  high 
momentum  regions  is  the  “outer”  velocity,  Uq,  we  can  use  the  equilibrium  logarithmic 
law [49]  and  aerodynamic  roughness  length,  z0,  to  predict  Uq\ 

Uo=  1 

Ur  K 

rmaY  (z)  is  the  advective  lag  between  passage  of  a  coherent  motion  at  elevation,  z,  and 
associated  modulation  of  processes  at  reference  height,  ZRef..  Thus,  Equations  9  and  10 


H 

zo 


(10) 


Figure  6:  (Color)  Sketch  of  low  momentum  re¬ 
gion  (LMR)  above  array  of  cubes.  The  LMR, 
u1  <  0,  is  denoted  by  transparent  gray.  Encap¬ 
sulating  the  LMR  are  hairpins  (blue).  The  LMR 
is  inclined  at  angle,  7,  and  the  hairpin  heads  ex¬ 
hibit  positive  transverse  vorticity  (sketch  denotes 
sweep  event,  Q4).  The  individual  hairpin  “legs” 
would  themselves  be  inclined  at  ~  45°  [41]. 
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can  be  combined  to  obtain  the  advective  lag: 


z  -  ZRef.  _  «(z  -  gRef.) 

tan  (7)  U0  tan  (7)  uT  log (H/ zq) 


(11) 


Normalizing  Tmax.(z)  by  friction  velocity  and  boundary  depth  (“shear  normalized”)  yields: 


Tmax.(z)uTH 


-1 


k(z-  2Ref.) 
tan  (7)  H  log  (if/ zq)  ' 


(12) 


We  have  used  the  LES  data  to  evaluate  z0,Eff.  for  the  block  cases  (h/H  =  1/8  and 
=  1/4),  enabling  prediction  of  the  outer  velocity,  Uq,  if  assuming  that  the  inertial  layer 
profile  retains  a  logarithmic  form  to  z  =  Zq.  Alternatively,  one  could  use  predictive 
models  for  z0  based  on  attributes  of  the  topography  [50,  51,  1]  although  this  would  require 
empirical  parameters [52],  For  Cases  SCI,  SC2,  and  Ul,  Figure  5  shows  predictions  from 
Equation  12  (solid  lines).  For  z^j  / H  <  z/H  <  0.2,  Equation  12  predictions  exhibit 
reasonable  agreement  with  the  LES  results  for  Case  Ul.  This  is  the  case  of  homogeneous 
roughness  in  which  the  roughness  sublayer  may  be  z/H  <  0.1.  For  Cases  SCI  and  SC2, 
Equation  12  predictions  agree  somewhat  well  with  the  LES  data  for  z^J H  <  z/H  < 
0.4,  corresponding  with  2  to  4  times  the  cube  height.  Thus,  the  model  performance  is 
moderately  successful  in  the  roughness  sublayer  and  into  the  inertial  layer. 


3  Dissemination 

The  work  described  above  (from  Year  #  1  of  the  research  effort)  has  been  presented 
through  invited  seminars,  conference  talks/proceedings,  and  peer  review  publications. 
Although  the  material  described  in  Section  1  (passive  scalar  fluxes  from  fractal  topog¬ 
raphy)  was  performed  before  the  funding  period  began,  we  include  dissemination  of  this 
material  below  as  part  of  the  “outcomes”  since  it  was  discussed  in  the  original  proposal. 

Invited  Talks : 

•  August,  2014:  Microscale  Atmospheric  Dynamics  Branch,  Army  Research  Labora¬ 
tory,  Adelphi,  MD. 

•  June,  2014:  Summer  Institute  on  Medicine  and  Energy,  Texas  Tech  University. 

•  March,  2014:  Mechanical  Science  and  Engineering  Department,  University  of  Illi¬ 
nois  at  Urbana-Champaign. 

•  March,  2014:  Mechanical  Engineering  Department,  University  of  Texas  at  Austin. 

•  March,  2014:  Sandia  National  Laboratory,  Albuquerque,  NM. 

•  February,  2014:  Mechanical  Engineering  Department,  University  of  Texas  at  San 
Antonio  [Declined]. 

•  January,  2014:  Mechanical  Engineering  Department,  University  of  Texas  at  Dallas. 

•  January,  2014:  Extreme  Fluids  Group,  Los  Alamos  National  Laboratory,  Los  Alamos, 
NM. 


•  July,  2013:  Army  Research  Laboratory,  Atmospheric  Dynamics  Branch. 

•  March,  2013:  Texas  Tech  University,  Department  Mechanical  Engineering,  Lub¬ 
bock,  TX. 

•  February,  2013:  Aerospace  Engineering  and  Engineering  Mechanics,  the  University 
of  Texas  at  Austin,  Austin,  TX. 

•  2012:  September,  2012:  Baylor  University,  Mechanical  Engineering  Department. 

Conference  Proceedings : 

•  Anderson  W,  Li  Q,  Bou-Zeid  E,  2014:  Proc.  of  American  Geophysical  Union,  Fall 
Meeting,  San  Francisco,  CA. 

•  Anderson  W,  Li  Q,  Bou-Zeid  E,  2014:  Proc.  of  American  Physical  Society,  Division 
of  Fluid  Dynamics,  San  Francisco,  CA. 

•  Li  Q,  Bou-Zeid  E,  Anderson  W,  2014:  Proc.  of  American  Physical  Society,  Division 
of  Fluid  Dynamics,  San  Francisco,  CA. 

•  Anderson  W,  Li  Q,  Bou-Zeid  E,  2014:  Proc.  of  American  Meteorological  Society, 
Symposium  on  Boundary  Layers  and  Turbulence,  Leeds,  England. 

•  Li  Q,  Bou-Zeid  E,  Anderson  W,  2014:  Proc.  of  American  Meteorological  Society, 
Atlanta,  GA. 

•  Anderson  W,  2012:  Proc.  of  American  Meteorological  Society,  Symposium  on 
Boundary  Layers  and  Turbulence,  Boston,  MA. 

•  Anderson  W,  Passalacqua  P,  2012:  Proc.  of  American  Meteorological  Society,  Sym¬ 
posium  on  Boundary  Layers  and  Turbulence,  Boston,  MA 

Peer-Review  Journal : 

•  Anderson  W,  2013:  Passive  scalar  roughness  lengths  for  atmospheric  boundary 
layer  flow  over  complex,  fractal  topographies.  Environmental  Fluid  Mechanics  13, 
479-501. 

Summary  :  This  manuscript  outlines  the  material  is  Section  1,  showing  that  Zqs/zq 
can  be  <<  10”1  for  multiscale,  fractal  topographies. 

•  Anderson  W,  Li  Q,  Bou-Zeid  E,  2014:  Transient  dynamics  of  coherent  motions  in 
an  urban-like  roughness  sublayer.  Journal  of  Turbulence  (Under  Review). 
Summary  :  This  manuscript  presents  the  material  briefly  described  in  Section  2. 
The  manuscript  was  submitted  on  July  10,  2014. 

4  Fiscal 

The  Year  ff  1  budget  was  $40,000,  which  supported  the  following  items: 

•  Travel  by  Anderson  to  the  American  Meteorological  Society,  Symposium  on  Bound¬ 
ary  Layers  and  Turbulence  in  Leeds,  UK 

•  Two  months  of  summer  salary  for  Anderson 
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•  Modest  equipment  and  publishing  budget 

In  Years  #  2  (FY15)  and  #  3  (FY16),  the  award  will  provide  $76084  and  $34840,  respec¬ 
tively.  These  funds  will  support  Anderson  summer  salary,  travel,  and  research  stipend 
for  a  new  doctoral  candidate  (Xiaowei  Zhu,  previously  a  research  associate  at  Nanyang 
Technological  University,  Singapore).  Comprehensive  details  of  expenses  and  fund  allo¬ 
cations  have  been  provided  on  the  full  proposal  submitted  to  ARO  in  2014,  as  part  of 
Anderson’s  new  start  at  UTD. 

5  Upcoming 

Year  #  2  of  the  project  (Year  #  1  of  the  UT  Dallas  new  start  award)  will  focus  on 
evaluating  the  transient  dynamics  of  flows  over  more  complex  topographies  (i.e.  beyond 
staggered,  uniform  height  cubes)  in  the  context  of  the  existing  Section  2  insights.  That 
is:  how  does  spacing  or  variable  height,  for  example,  affect  the  Equation  9  to  12  de¬ 
velopment.  In  fact,  the  topic  of  spacing  is  already  “underway”  through  a  collaboration 
between  Anderson  and  Elic  Bou-Zeid,  Princeton  University.  Elie  and  his  student,  Qi 
Li,  have  studied  the  transient  dynamics  questions  with  spacing  varied  above  and  below 
the  cases  considered  here.  In  addition,  we  will  begin  to  consider  the  effects  of  fractal 
topographies  (with  varying  spectral  exponent,  or  fractal  dimension),  to  evaluate  how 
varying  “roughness”  modifies  the  spatial  extent  of  coherent  parcels  of  fluid.  Townsend’s 
outer  layer  similarity  hypothesis  states  that  smooth  and  rough  wall  turbulence  statistics 
should  be  invariant  in  the  outer  layer  [53].  Some  recent  efforts  have  demonstrated  that 
complex  roughness  serves  to  attenuate  coherence  in  the  roughness  sublayer [5,  42],  while 
Hong  et  al.  [54]  demonstrated  a  persistent  “near  wall”  signature  in  outer  layer  turbu¬ 
lence  statistics  for  flow  over  a  distribution  of  pyramids.  Thus,  by  considering  a  suite  of 
fractal  topographies  with  varying  fractal  dimension,  we  could  determine  the  presence  of 
monotonic  trends  in  turbulence  morphology  with  additional  geometric  complexity  of  the 
underlying  topography. 

6  Institutional  Move:  Baylor  University  to  the  University  of  Texas  at  Dallas 

Anderson  completed  his  PhD  (Mechanical  Engineering)  at  The  Johns  Hopkins  University 
in  July,  2011,  and  began  as  a  tenure-track  faculty  in  the  Mechanical  Engineering  Depart¬ 
ment  at  Baylor  University  in  Fall,  2011.  Anderson  greatly  enjoyed  being  at  Baylor,  but 
experienced  acute  difficulty  recruiting  graduate  students.  This  was  exacerbated  by  virtue 
of  federal  grants  awarded  and  a  growing  research  portfolio  (please  see  Section  7);  in  fact, 
Anderson  in  one  case  had  to  decline  significant  funds  because  they  could  never  have  been 
expended.  Thus,  Anderson  was  in  the  unfortunate  situation  of  having  to  move  to  an 
institution  with  a  more  matured  research  enterprise.  In  July,  2014,  Anderson  moved  to 
the  Mechanical  Engineering  Department  at  The  University  of  Texas  at  Dallas. 

7  Additional  Funding 

In  the  interest  of  openness,  Anderson  has  below  provided  a  list  of  active  and  pending 
grants,  and  a  brief  summary  of  the  research.  Anderson  will  also  gladly  provide  contact 
information  for  the  respective  program  manager  at  ARO’s  request'. 

'Note  that  this  list  excludes  white  papers  presently  under  review  or  external  consulting  activities. 
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•  $240000:  Air  Force  Office  of  Scientific  Research,  Young  Investigator  Program  (Oc¬ 
tober,  2014  -  September,  2017) 

Summary.  This  AFOSR  award  will  support  a  new  topic  studying  the  role  of 
canonical  roughness  in  sustaining  mean  secondary  flows  within  a  turbulent  bound¬ 
ary  layer. 

•  $20000:  Shell  Exploration  and  Production  Company  (January,  2014  -  December, 
2014) 

Summary.  This  industrial  grant  supports  efforts  to  model  atmospheric  surface 
layer  flows  over  fields  of  aeolian  sand  dunes.  The  work  is  motivated  by  the  need 
to  understand  aerodynamic  stresses  responsible  for  geomorphic  evolution  of  dune 
fields. 

•  $1.2  million:  Air  Force  Office  of  Scientific  Research,  Turbulence  and  Transition 
Program  (July,  2014  -  June,  2018) 

Summary.  This  award  is  to  be  divided  between  Anderson  and  collaborators  (K. 
Christensen,  Notre  Dame,  and  C.  Pantano,  UIUC).  The  award  focuses  on  refine¬ 
ment  of  wall  models  for  large-eddy  simulation  of  rough  wall  turbulent  boundary 
layers.  The  refinement  here  will  be  accomplished  by  adaptation  of  ‘Amplitude 
Modulation’  concepts  -  developed  recently  for  smooth  wall  flows  by  Ivan  Marusic 
and  collaborators  -  for  flows  over  rough  surfaces. 

•  Under  Review:  National  Science  Foundation,  CAREER  (Fluid  Dynamics  Program, 
CBET)  (September,  2015-August,  2020) 

Summary  :  This  project  will  support  a  more  detailed  study  of  turbulence  processes 
in  dune  field  environments  (that  is,  development  of  bulk  characterizations  of  the 
turbulent  mixing- layer  like  processes  present  in  close  proximity  to  the  dunes).  This 
work  will  especially  focus  on  “realistic”  aeolian  dune  fields,  with  spatial  distribution 
set  by  winds  from  multiple  directions. 

•  Under  Review:  National  Science  Foundation,  Physical  and  Dynamic  Meteorology 
Program  (GEO) 

Summary:  Anderson  has  submitted  a  DoD-style  white  paper  to  this  program 
for  consideration.  At  the  program  manager’s  request,  Anderson  is  now  developing 
a  full  proposal.  The  project  will  focus  on  the  morphology  of  coherent  turbulent 
parcels  in  the  convectively  stratified  atmospheric  boundary  layer,  and  their  role  in 
intermittent  surface  fluxes  of  dust  from  arid  landscapes  on  the  Llano  Estacado  in 
west  Texas  and  eastern  New  Mexico. 

•  Linder  Review:  National  Science  Foundation,  Fluid  Dynamics  Program  (CBET) 
Summary:  Anderson  (as  co-PI)  and  collaborators  (PI  K.  Christensen,  University 
of  Notre  Dame;  co-PI  Gianluca  Blois,  University  of  Notre  Dame;  co-PI  J.  Best, 
University  of  Illinois  at  Urbana-Champaign;  co-PI  G.  Kocurek,  University  of  Texas 
at  Austin),  have  submitted  a  full  proposal  to  the  NSF-ENF-CBET-Fluid  Dynamics 
Program  that  focuses  on  turbulent  flows  over  canonical  barchan  dunes.  The  pro¬ 
posal  is  an  experimental-numerical  collaborative.  The  three-year  budget  request 
was  $540,000,  of  which  Anderson  will  receive  $180,000. 

•  Linder  Review:  Texas  Department  of  Transportation  (TXDoT),  Research  and  Tech¬ 
nology  Implementation  Office  (RTI) 
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Summary.  A  three-year  project  is  under  review  at  TXDoT  with  Anderson  as  the 
sole  PI.  The  project  will  use  the  Weather  Research  and  Forecasting  (WRF)  com¬ 
putation  package  for  modeling  mesoscale  meteorological  conditions  and  using  short 
term  forecasts  to  inform  trafficability  in  Texas. 

•  To  be  Submitted:  Office  of  Naval  Research,  Young  Investigator  Program  (Program¬ 
matic  Area:  Hull  Performance/Undersea  Hydromechanics) 

Summary.  This  research  is  complimentary  to  aspects  of  this  ongoing  ARO  project. 
In  this  effort,  the  LES  channel  flow  code  will  be  used  to  model  flow  over  a  suite 
of  complex,  fractal  (and  multifractal)  topographies  to  evaluate  how  modifying  geo¬ 
metric  attributes  of  the  synthetic  topographies  affects  the  turbulence  statistics. 

8  Army  Research  Laboratory  Collaborations 

Anderson  visited  the  Army  Research  Laboratory  Headquarters  in  Adelphi,  MD  in  July, 
2013  and  August,  2014.  The  2014  visit  was  also  intended  as  a  “status  update”  to  staff 
scientists  there  who  were  keenly  aware  of  Anderson’s  research  (Drs.  Yansen  Wang,  Ben 
MacCall,  Cheryl  Klipp,  Chatt  Williamson,  and  Sandra  Collier).  In  addition  to  discussing 
Anderson’s  activities,  there  is  always  in-depth  discussion  about  non-secure  research  being 
conducted  there.  Anderson  and  Ben  MacCall  continue  to  explore  potential  collaboration 
opportunities  in  diverse  turbulence  problems. 

References  Cited 

[1]  W.  Anderson  and  C.  Meneveau.  A  dynamic  large-eddy  simulation  model  for  bound¬ 
ary  layer  flow  over  multiscale,  fractal-like  surfaces.  J.  Fluid  Mech.,  679:288-314, 
2011. 

[2]  W.  Anderson,  P.  Passalacqua,  F.  Porte-Agel,  and  C.  Meneveau.  Large-eddy  simula¬ 
tion  of  atmospheric  boundary  layer  flow  over  fluvial-like  landscapes  using  a  dynamic 
roughness  model.  Boundary- Layer  Meteorol.,  144:263-286,  2012. 

[3]  I.  Rodrfguez-Iturbe  and  A.  Rinaldo.  Fractal  river  basins:  chance  and  self¬ 
organization.  Cambridge  University  Press,  Cambridge,  USA,  1997. 

[4]  Z.-T.  Xie,  O.  Coceal,  and  I.P.  Castro.  Large-eddy  simulation  of  flows  over  random 
urban-like  obstacles.  Boundary- Layer  Meteorol .,  129:1-23,  2008. 

[5]  O.  Coceal,  A.  Dobre,  T.  G.  Thomas,  and  S.E.  Belcher.  Structure  of  turbulent  flow 
over  regular  arrays  of  cubical  roughness.  J.  Fluid  Mech.,  589:375-409,  2007. 

[6]  Margi  Bohm,  John  J.  Finnigan,  Michael  R.  Raupach,  and  Dale  Hughs.  Turbulence 
structure  within  and  above  a  canopy  of  bluff  elements.  Boundary  Layer  Meteorol., 
146:393-419,  2013. 

[7]  M.  Bohm,  J.  J.  Finnigan,  and  M.R.  Raupach.  Dispersive  fluxes  and  caopy  flows:  Just 
how  important  are  they?  Proc.  2fth  Conf.  on  Agr.  and  Forest  Meteorol.,  American 
Meteorol.  Soc.,  Davis,  Calif.,  pages  106-107,  2000. 

[8]  M.R.  Raupach,  R.A.  Antonia,  and  S.  Rajagopalan.  Rough-wall  turbulent  boundary 
layers.  Appl.  Mech.  Rev.,  44:1-25,  1991. 

[9]  S.E.  Belcher,  I.N.  Harman,  and  J.J.  Finnigan.  The  wind  in  the  willows:  flows  in 
forest  canopies  in  complex  terrain.  Ann.  Rev.  Fluid  Mech.,  44:479-504,  2012. 

[10]  J.  Finnigan.  Turbulence  in  plant  canopies.  Ann.  Rev.  Fluid  Mech.,  32:519-571,  2000. 


12 


[11]  M.  Germano,  U.  Piomclli,  P.  Moin,  and  W.  Cabot.  A  dynamic  subgrid-scale  eddy 
viscosity  model.  Phys.  Fluids  A,  3:1760-1765,  1991. 

[12]  A.S.  Monin  and  A.M.  Obukhov.  Basic  laws  of  turbulent  mixing  in  the  ground  layer 
of  the  atmosphere.  Tr.  Geofiz.  Inst.,  Akad  Nauk  SSSR,  151:163-187,  1954. 

[13]  Vijayant  Kumar,  Jan  Klcissl,  Charles  Meneveau,  and  Marc  B.  Parlange.  Large- 
eddy  simulation  of  a  diurnal  cycle  of  the  atmospheric  boundary  layer:  Atmospheric 
stability  and  scaling  issues.  Water  Resour.  Res.,  42:W06D09,  2006. 

[14]  H.U.  Sverdrup.  On  the  evaporation  from  the  oceans.  J.  Marine  Res.,  1:3-14,  1937. 

[15]  W.  Brutsaert.  A  model  for  evaporation  as  a  molecular  diffusion  process  into  a 
turbulent  atmosphere.  J.  Geophys.  Res.,  70:5017-5024,  1965. 

[16]  W.  Brutsaert.  A  theory  for  local  evaporation  (or  heat  transfer)  from  rough  and 
smooth  surfaces  at  ground  level.  Water  Resour.  Res.,  11:543-550,  1975. 

[17]  P.R.  Owen  and  W.R.  Thomson.  Heat  transfer  across  rough  surfaces.  J.  Fluid  Mech., 
15:321-334,  1963. 

[18]  N.  Sheriff  and  P.  Gumlcy.  Heat-transfer  and  friction  properties  of  surfaces  with 
discrete  roughness.  Int.  J.  Heat  and  Mass  Transfer,  9:1297-1320,  1966. 

[19]  D.F.  Dipprey  and  R.M.  Sabersky.  Heat  and  momentum  transfer  in  smooth  and  rough 
tubes  at  various  prandtl  numbers.  Int.  J.  Heat  Mass  Transfer ,  6:329-332,  1963. 

[20]  A.C.M.  Beljaars  and  A. A.M.  Holtslag.  Flux  parameterization  over  land  surfaces  for 
atmospheric  models.  J.  Appl.  Meteor.,  30:327-341,  1991. 

[21]  A.T.  Cahill,  M.B.  Parlange,  and  J.D.  Albertson.  On  the  brutsaert  temperature 
roughness  length  model  for  sensible  heat  flux  estimation.  Water  Resour.  Res., 
33:2315-2324,  1997. 

[22]  A.M.  Yaglom  and  B.A.  Kader.  Heat  and  mass  transfer  between  a  rough  wall  and 
turbulent  fluid  flow  at  high  reynolds  and  peclet  numbers.  J.  Fluid  Mech.,  62:601-623, 
1974. 

[23]  W.  Brutsaert.  Evaporation  into  the  atmosphere.  Kluwer  Academic  Publishers,  Nor- 
well,  MA,  1982. 

[24]  W.  Anderson.  Passive  scalar  roughness  lengths  for  atmospheric  boundary  layer  flow 
over  complex,  fractal  topographies.  Env.  Fluid  Mech.,  pages  DIO  10.1007/sl0652- 
013-9272-9,  2013. 

[25]  W.  Anderson.  An  immersed  boundary  method  wall  model  for  high-reynolds  number 
channel  flow  over  complex  topography.  Int.  J.  Numer.  Methods  Fluids,  71:1588-1608, 
2012. 

[26]  M.  Calaf,  M.B.  Parlange,  and  C.  Meneveau.  Large  eddy  simulation  study  of  scalar 
transport  in  fully  developed  wind-turbine  array  boundary  layers.  Phys.  Fluids, 
23:126603-16,  2011. 

[27]  J.J.  Finnigan,  R.H.  Shaw,  and  E.G.  Patton.  Turbulence  structure  above  a  vegetation 
canopy.  J.  Fluid  Mech.,  637:387-424,  2009. 

[28]  H.  Cheng  and  I.P.  Castro.  Near  wall  flow  over  urban-like  roughness.  Boundary- Layer 
Meteorol,  104:229-259,  2002. 

[29]  S.E.  Belcher,  N.  Jerram,  and  J.C.R.  Hunt.  Adjustment  of  a  turbulent  boundary 
layer  to  a  canopy  of  roughness  elements.  J.  of  Fluid  Mech.,  488:369-398,  2003. 

[30]  W.  Anderson  and  M.  Chamecki.  Physical  Review  E,  89:013005-1-14,  2014. 

[31]  Ian  Harmon  and  John  J.  Finnigan.  A  simple  unified  theory  for  flow  in  the  canopy 
and  roughness  sublayer.  Boundary  Layer  Meteorol. ,  123:339-364,  2007. 

[32]  I.P.  Castro.  Rough- wall  boundary  layers:  mean  flow  universality.  J.  Fluid  Mech., 
585:469-485,  2007. 


13 


[33]  M.R.  Raupach,  J.J.  Finnigan,  and  Y.  Brunet.  Coherent  eddies  and  turbulence  in 
vegetation  canopies:  the  mixing  layer  analogy.  Boundary- Layer  Meteorol,  78:351- 
382,  1996. 

[34]  M.  Ghisalberti.  Obstructed  shear  flows:  similarities  across  systems  and  scales.  J. 
Fluid  Meek,  641:51,  2009. 

[35]  M.M.  Rogers  and  R.D.  Moser.  Direct  simulation  of  a  self-similar  turbulent  mixing 
layer.  Phys.  Fluids  A,  6:903-922,  1994. 

[36]  N.  Hutchins  and  I.  Marusic.  J.  Fluid  Mech .,  579:1-28,  2007. 

[37]  J.  Zhou,  R.  J.  Adrian,  S.  Balachandar,  and  T.M.  Kendall.  Mechanisms  for  generating 
coherent  packets  of  hairpin  vortices  in  channel  flow.  J.  Fluid  Mech.,  387:353-359, 
1999. 

[38]  R.J.  Adrian,  C.D.  Meinhart,  and  C.D.  Tomkins.  Vortex  organization  in  the  outer 
region  of  the  turbulent  boundary  layer.  J.  Fluid  Mech.,  422:1-54,  2000. 

[39]  K.T.  Christensen  and  R.J.  Adrian.  Statistical  evidence  of  hairpin  vortex  packets  in 
wall  turbulence.  J.  Fluid  Mech.,  431:433-443,  2001. 

[40]  B.  Ganapathisubramani,  E.  K.  Longmire,  and  I.  Marusic.  J.  Fluid  Mech.,  478:35-46, 
2003. 

[41]  R.J.  Adrian.  Phys.  Fluids,  19:041301,  2007. 

[42]  Y.  Wu  and  K.  T.  Christensen.  Spatial  structure  of  a  turbulent  boundary  layer  with 
irregular  surface  roughness.  J.  Fluid  Mech.,  655:380-418,  2010. 

[43]  D.J.C.  Dennis  and  T.B.  Nickels.  J.  Fluid  Mech.,  673:180-217,  2011. 

[44]  D.J.C.  Dennis  and  T.B.  Nickels.  J.  Fluid  Mech.,  673:218-244,  2011. 

[45]  C.D.  Meinhart  and  R.J.  Adrian.  On  the  existence  of  uniform  momentum  zones  in  a 
turbulent  boundary  layer.  Phys.  Fluids,  7:694-696,  1995. 

[46]  R.J.  Adrian.  Linking  correlations  and  structure:  stochastic  estimation  and  condi¬ 
tional  averaging.  In  Zoran  P.  Zaric  Memorial  International  Seminar  on  Near-  Wall 
Turbulence.  Hemisphere,  1988. 

[47]  W.  Anderson,  J.  Barros,  and  K.  T.  Christensen.  J.  Fluid  Mech.,  page  Linder  Rev., 
2014. 

[48]  S.S.  Lu  and  W.W.  Willmarth.  Measurements  of  the  structure  of  reynolds  stress  in  a 
turbulent  boundary  layer.  J.  Fluid  Mech.,  60:481-571,  1973. 

[49]  Lb  Piomelli  and  E.  Balaras.  Wall-layer  models  for  large-eddy  simulation.  Ann.  Rev. 
Fluid  Mech. ,  34:349-374,  2002. 

[50]  J.R.  Garratt.  The  atmospheric  boundary  layer.  Cambridge  LIniversity  Press,  LIK, 
1994. 

[51]  M.V.  Zagarola  and  A.J.  Smits.  Mean-flow  scaling  of  turbulent  pipe  flow.  J.  Fluid 
Meek,  373:33-79,  1998. 

[52]  M.P.  Schultz  and  K.A.  Flack.  Turbulent  boundary  layers  on  a  systematically-varied 
rough  wall.  Phys.  Fluids,  21:015104,  2009. 

[53]  A.  A.  Townsend.  The  structure  of  turbulent  shear  flow.  Cambridge  LIniversity  Press, 
Cambridge,  UK,  1976. 

[54]  J.  Hong,  J.  Katz,  C.  Meneveau,  and  M.  Schultz.  Coherent  structures  and  associated 
subgrid-scale  energy  transfer  in  a  rough-wall  channel  flow.  J.  Fluid  Mech.,  712:92- 
128,  2012. 


14 


n 


4 

2 

0 

-2 

-4 


n 


2 

0 

-2 


(c)  1 

^0.5 


n 


■  it 8  Q 


4 


0 


0.05 


0.1 


0.15 

tuTH~l 


~02 


0.25 


(d)  1 

i£! 


0.5-  , 


,.v  a 


K 


> 


0.05 


_ i _ 

0.1 


.  ©' 


0.15 

tuTH~l 


0.2 


i 

0.25 


~03 


n 


(<0 

tci 


0.5 


I 


11 


«  «. , ' 


5 

0 

U-5 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8 

tuTH~l 


(f) 


3  0.5P 


I 


111 


01 


0 


U-2 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8 

tuTH~l 


Figure  4:  (Color)  Time-height  contours  of  fluctuating  velocity  (u'(x,t)  =  u(x,t)  — 
(u(x,t))t)  components  for  flow  over  cubical  topography  (Panels  a-d)  and  homogeneous 
roughness  (Panels  e  and  f)  at  positions  indicated  in  Figure  3:  (a)  u'(xi,yi,z,t)  at  x\\  (b) 
w'(xi,yi,  z,t)  at  x\\  (c)  u'(x^y^z,t)  at  x4;  (d)  w'(x^y^,z,t)  at  x 4;  (e)  u'(xc,yc,z,t )  at 
center;  and  (f)  w'(xc,yc,z,t)  at  center.  In  addition  to  color  floods,  line  contours  denote 
contribution  to  u'w'  due  to  Q 4  “sweep”  (black)  and  Q 2  “ejection”  (yellow)  events.  Note 
that  in  Parts  (a)  and  (b),  u'(x,y,z/ H  <  1/4,  t)  =  0  owing  to  the  presence  of  cubes  of 
height  h/H  =  1/4 
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Figure  5:  (Color)  Advective  lag  versus  wall-normal  elevation  for  Cases  U1  (blue  symbols), 
SCI  (black  symbols),  and  SC2  (red  symbols).  For  Cases  SCI  and  SC2,  symbols  corre¬ 
spond  with:  xi  (circles),  x2  (squares),  x3  (“plus”  sign),  and  x4  (asterisk).  Blue  symbols 
correspond  with  xc  for  the  homogeneous  rough  case.  Vertical  dashed  lines  denote  ZRef.  for 
different  cases  (dashed  blue:  Ul;  dashed  black:  SCI;  and  dashed  red:  SC2).  Solid  lines 
represent  Equation  12  predictions  of  time-lag  for  effective  roughness  lengths  determined 
a  posteriori  from  simulations. 
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